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Abstract 

Activated surface diffusion with interacting adsorbates is analyzed within the Linear Response 
Theory framework. The so-called interacting single adsorbate model is justified by means of a 
two-bath model, where one harmonic bath takes into account the interaction with the surface 
phonons, while the other one describes the surface coverage, this leading to defining a collisional 
friction. Here, the corresponding theory is applied to simple systems, such as diffusion on flat 
surfaces and the frustrated translational motion in a harmonic potential. Classical and quantum 
closed formulas are obtained. Furthermore, a more realistic problem, such as atomic Na diffusion 
on the corrugated Cu(OOl) surface, is presented and discussed within the classical context as well 
as within the framework of Kramer's theory. Quantum corrections to the classical results are also 
analyzed and discussed. 
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I. INTRODUCTION 



The main purpose of spectroscopic experiments involving a probe and a system at thermal 
equilibrium with a reservoir or thermal bath consists of measuring the system response 
under the perturbation caused by the probe . The intrinsic properties of matter can 
be extracted by analyzing this response, which thus becomes a very important tool for 
physicists, who can elucidate the microscopic structure and its dynamics through it. Many 
times this response can be described by first order perturbation theory and determined by 
the spectrum of the spontaneous fluctuations of the reservoir, as established by the so-called 
fluctuation-dissipation (FD) theorem jsl- Within this linear approximation, the probability 
per time unit that the full system formed by probe and reservoir changes from the initial state 
to the final state is given by the Fermi golden rule. More specifically, the probe transition 
probability from an initial state to a certain final state is given by the time Fourier transform 
of the autocorrelation function associated with the operator defined by the matrix element 
of the corresponding interaction Hamiltonian. This type of studies is based on the work 
developed by van Hove 0, [sl, who introduced the space-time correlation function (namely 
the van Hove function), a generalization of the well-known pair-distribution function from 
the theory of liquids, as a tool to study the scattering of probe particles (slow neutrons) by 
quantum systems consisting of interacting particles at thermal equilibrium. Within the Born 
approximation in scattering theory, the nature of the scattered particles as well as the details 
of the system-probe interaction potential are largely irrelevant, this essentially reducing 
the scattering problem to a typical statistical mechanics problem {g]. The linear response 
function of a system consisting of interacting particles, also known as dynamic structure 
factor or scattering law, can be then related to the spontaneous-fluctuation spectrum of 
such particles and expressed in terms of particle density-density correlation functions js-S]. 
As is known, a complete description of this dynamics can also be obtained through the 
linear response theory, where the FD theorem is used to derive alternative response function 
expressions for the dynamic structure factor. 

From quasielastic He atom scattering (QHAS) experiments at low energies very detailed 
information about defect and adsorbate dynamics on surfaces can be obtained IJQ. Man- 



son and Celli 



10| generalized van Hove's theory of neutron scattering by crystals and liquids 



to atom surface scattering within the transition matrix formalism. Within this approach. 



2 



small coverage of adsorbates and/or defects were assumed in order to ignore both their 
interactions as well as multiple scattering with He atoms. Diffuse elastic and inelastic fea- 
tures were interpreted again through the dynamical structure factor which is proportional 
to the observed line shapes and is expressed in terms of the transfer of energy and parallel 
momentum to the He atoms before and after the scattering process. The QHAS technique 
has also been applied to study surface diffusion on metals with different types of atomic 
and molecular adsorbates, the diffusion of Na adatoms (at different coverages) on Cu(OOl) 
being one of the most extensively studied systems. In the case of massive particles, where 
very large timescales are involved, the dynamic structure factor can be expressed in terms of 
the adsorbate positions, which has led to the so-called single adsorbate model. Within this 
approach, where low coverage is assumed, surface diffusion (i.e., the adsorbate motion) is 
described by classical stochastic trajectories of adsorbates issued from solving the standard 
Langevin equation. As in the case of Brownian-like particles, this equation encompasses two 
contributions: (1) a (external) force arising from the deterministic, phenomenological adia- 
batic potential describing the adsorbate-substrate interaction at zero surface temperature; 
and (2) a stochastic force (usually, a Gaussian white noise) accounting for the vibrational 
effects induced by the temperature on the surface lattice atoms (and therefore on the adsor- 



bates). The system dynamics is then obtained after solving the Langevin equation 
ana^ 
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yzing the results derived from it in terms of the so-called motional narrowing effect 
16] as well as Kramers' turnover theory 17, 18\ and the dephasing theory [3]. Within 



this scenario, it is usually assumed that He-adsorbate interactions play no role on the surface 
dynamics, thus being scarcely analyzed. 

When the coverage is increased, the dynamic structure factor also provides valuable 
information about the nature of the adsorbate-adsorbate interaction, which should be in- 
cluded in the corresponding theoretical studies. In this way, pairwise potential functions 
accounting for the adsorbate-adsorbate interactions are usually introduced into Langevin 

n 

molecular djTiamics (LMD) simulations ^J]. Recently, an alternative procedure has been 
considered, where such pairwise interactions are described by a purely stochastic model, 
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423l|. This approach, also within 



namely the interacting single adsorbate (ISA) model 
the standard Langevin framework, is based on both the theory of spectral-line collisional 
broadening developed by van Vleck and Weisskopf 2J] and the elementary kinetic theory 
of gases 2(], and explains fairly well the experimental broadening observed with increasing 
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coverage [l^. The standard Langevin equation is solved with two different, non-correlated 
noise functions: (1) a Gaussian white noise accounting for the surface friction, as before, and 
(2) a white shot noise |25| replacing the pairwise interaction potential which simulates the 
adsorbate-adsorbate collisions. A double Markovian assumption is therefore considered (for 
the interaction with the surface and for the interaction among adsorbates). This assumption 
holds because, on the one hand, substrate excitation timescales are much shorter than the 
timescales associated with the adatom motions (the maximum frequency of the substrate 
excitation is around 20-30 meV, while the characteristic vibrational frequency of the adatom 
is around 4-6 meV). On the other hand, the time between consecutive collisions (measured 
through a collisional frequency, which would depend on temperature, coverage or the adpar- 
ticle mass) is typically much longer than the effective time that a collisional effect may last 
(i.e., the time two adsorbates may effectively be in physical contact). Therefore, memory 
effects are not taken into account. Just as the adsorbate-surface interaction is characterized 
by a (surface) friction, the adsorbate-adsorbate interaction will also be describable in terms 
of a collisional friction, which will vary with coverage. With this simple stochastic model, 
where the total friction is the sum of the substrate friction and the collisional friction, a 
good agreement with the experimental results for coverages up to 0.12, approximately, has 
been obtained. For higher coverage values, the model cannot be applied due to the appear- 
ance of ordered structures, as it has been observed [26] experimentally between 0.12 and 
0.16. Recently, the collisional friction has been estimated from experiments with benzene 
on graphite |27j. Although further investigation at a microscopic level and first-principle 
calculations are needed, this simple stochastic model at moderate coverages is able to pro- 
vide a complementary view of diffusion (through the quasielastic Q-peak) and low frequency 
vibrational motions (through the frustrated translational T-mode peak), at or around zero 
energy transfers (very long time dynamical processes), respectively. This could be under- 
stood because any trace of the true interaction potential seems to be wiped out due to the 
relatively large number of collisions taking place at very long times (the time scale for the 
diffusion regime). Actually, this purely stochastic model can be derived from a microscopic 
classical Hamiltonian model characterized by two baths associated with two independent 
collections of harmonic oscillators 28| , one related to the surface phonons and the other one 
describing the presence of adsorbates, the source of the collisional friction. In this model, 
the coupling to low-lying electron-hole excitations (electroninc friction) was not considered 
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for simplicity. One of the purposes in the present work is to extend this two-bath model to 
quantum activated diffusion. 

Surface diffusion processes are also studied experimentally by means of the so-called spin 
echo techniques. Thus, since the advent of He spin echo (HeSE) spectroscopy 29N31| and the 
improved signal in the neutron spin echo (NSE) spectrometer 32|, |33| , fast diffusion processes 
are more accessible and allow us to to better determine interaction potentials. In both types 
of experiments, the observable is the so-called intermediate scattering function or polarization 



function, which is the inverse time Fourier transform of the dynamic structure factor. 



nis IS 



a complex function whose real and imaginary parts can be observed experimentally |30l. 131 1. 
The intermediate scattering function is also the space Fourier transform of the van Hove 
function, which, in general, is also a complex-valued function. The complex character of 
these functions can be understood as a signature of the quantum nature of the diffusion 
dynamics. In this regard, a quantum Markovian theory of surface diffusion for interacting 
adsorbates has also been proposed recently [3^. The imaginary part of the van Hove function 
is important at small values of time or high temperatures (timescales of the order of hP, 
with /3 = l/ksT, ks being Boltzmann's constant). This dynamical regime takes place when 



the mean de Broglie wavelength = h/\/2mkBT (m is the adsorbate mass) is of the order 
of or greater than typical interparticle distances. For these timescales and distances, the 
adparticle positions are no longer variables, but Heisenberg operators that do not commute 
at two different times. This theory thus tries to reconcile classical and quantum calculations 
when no diffusion by tunneling is considered. 

The purpose of this work is to provide an analysis of activated surface diffusion with inter- 
acting adsorbates within the linear response theory framework. The theoretical framework 
of the aforementioned two-bath model is then applied to simple systems, such as diffusion 
on fiat surfaces and the frustrated translational motion in a harmonic potential, obtaining 
classical and quantum closed formulas. Moreover, a more realistic problem, such as atomic 
Na diffusion on on the corrugated Cu(OOl) surface, is presented and discussed within the 
classical context as well as within the framework of Kramer's theory, introducing quan- 
tum corrections to the classical results which will be analyzed and discussed. Note that 
an appropriate understanding of this surface dynamics is very important, for diffusion is a 
preliminary step in more complicated surface phenomena, such as heterogeneous catalysis, 
crystal growth, lubrication, associative desorption, etc. Furthermore, the QHAS and HeSE 
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techniques can be considered as the surface science analogue of the quasielastic neutron 
scattering techniques, which has been widely and successfully applied to analyze diffusion 
in bulk. 

According to our purposes, we have organized this work as follows. In Section [ITl we give 
a general overview of the Linear Response Theory applied to activated surface diffusion. In 
particular, a revision of the classical two-bath model is presented. Regarding the quantum 
version of this model, a proposal and discussion are also given in order to describe the 
stochastic trajectories of the adsorbates at different coverages. In Section IIIIl applications 
to simple models (fiat surfaces and driven damped harmonic oscillator) as well as to more 
realistic problems, as atomic Na diffusion on Cu(OOl), are analyzed within the context of 
classical and quantum dynamics. Kramers's turnover theory is also analyzed within the 
classical context. Finally, in Section IIVI we summarize the conclusions derived from this 
work as well as some future work. 



II. GENERAL THEORY FOR ACTIVATED SURFACE DIFFUSION 
A. Dynamic structure factor and intermediate scattering function 



Space-time correlation functions [35|] can be used to describe the decay of spontaneous 
thermal fluctuations at surfaces, being central to the study of transport phenomena. These 
functions are defined as the thermodynamic average of the product of two dynamical vari- 
ables, each one expressing the instantaneous deviation from its corresponding equilibrium 
value at particular points on the surface and time. A complete description of the particle 
dynamics in a many-body system is then reached when the behavior of the corresponding 
correlation functions over the entire wavenumber range is studied. This range splits into 
different characteristic regions, each one associated with a different set of properties of the 
system. In the case of scattering experiments, since the momentum and energy transfers are 
the relevant quantities, any correlation-function theory has to be developed necessarily in 
terms of such quantities. Space-time correlation functions can also be used to describe the 
response of a fluid under a weak, external perturbation. Indeed, the reason why space-time 
correlation functions are central quantities in transport phenomena in fluids is, precisely, 
because of the equivalence between spontaneous fluctuation and linear response. 
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In general, a surface local dynamical variable is defined as 

1 ^ 

A(R, t) = -= V a,(t)5(R - R,(t)), (1) 

.=1 

where ai(t) is any physical quantity and Ri(t) = {xi(t) , yi{t)) is the time-dependent center 
of the position operator of the adparticle on a two-dimensional surface. The corresponding 
fluctuation is usually defined as 6A(Tl,t) = A(Tl,t) — (A(R, t))/3, where the average on a 
canonical ensemble is denoted by ( ■ )^. The dynamical variable conserves if it satisfies a 
continuity equation of the form 

i(R,t) = -VR-JA(R,t), (2) 

where 3a is the current associated with the A variable and the dot over A denotes the total 
time derivative. Here, the dynamical variables of particular interest are the number density, 

1 ^ 

p(R,t) = -=5^5(R-R,(t)), (3) 

^ i=l 

and the current density 

1 ^ 

J(R,t) = ^v,,(t)5(R-R,(t)), (4) 



1=1 

where Vj accounts for the velocities of the N adparticles. The corresponding van Hove 
fluctuation density autocorrelation function 4, |6| reads as 

G'(|R-R'|,t) = S(5p(R',0)5p(R,t))/3, (5) 

and a similar expression holds for the current density. The adparticle density is given by 
p = N/T,, where S is the surface area and the coverage is defined by 6* = A^/A^max, with A^max 
being the maximum number of sites in the S area. In most physical systems, correlation 
effects are negligible at large space or time separation, the asymptotic limit being a simple 
product of thermodynamically averaged quantities. 



the observable 



In analogy to scattering of slow neutrons by crystals and liquids 
magnitude in QHAS experiments is the so-called differential reflection coefficient, 

-^^^ = n,7S{A^M. (6) 
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This coefficient gives the probabihty that the He atoms (probe particles) scattered from the 
interacting adsorbates on the surface reach a certain sohd angle Q with an energy exchange 
hu = Ef — Ei and wave vector transfer parallel to the surface AK = Kj — Kj. In Eq. (E]), 
rifi is the concentration of adparticles; J-" is the atomic form factor, which depends on the 
interaction potential between the probe atoms in the beam and the adparticles on the surface; 
and S'(AK, u) is the dynamic structure factor, which gives, apart from other peaks, the Q 
and T-mode peaks, also providing a complete information about the dynamics and structure 
of the adsorbates through particle distribution functions. Experimental information about 
long distance correlations is obtained from the dynamic structure factor when considering 
small values of AK, while information on long time correlations is provided at small energy 
transfers, hw. 

The pair distribution functions are given by means of the van Hove or time-dependent 
pair correlation function G{Ti, t) [4]. This function is related to the dynamic structure factor 
by a double Fourier transform, in space and time, as 

5(AK,a;) = ^ j j G(R, t)e^(^^-^-^*) rfRrft. (7) 

Given an adparticle at the origin at some arbitrary initial time, (^(R, t) represents the 
average probability to find a particle (the same or another one) at the surface position 
R = (x, y) at a time t. Thus, this function generalizes the well-known pair distribution 
function gijV) from statistical mechanics 2|,|7| by providing information about the interacting 
particle dynamics. 

The position operators of the adsorbates are given, in general, by the corresponding 
Heisenberg operators (defined for all j = 1, ■ ■ ■ , iV adparticles and time t), 

R.(t) = e'^*/''R,■e-'^*/^ (8) 

where H is the Hamiltonian of the total system. As mentioned above, the space Fourier 
transform of the G-function is the intermediate scattering function, 

/(AK,t) = N [[ G{R,t)e'^^-^dR 



^(pAK(t)piK(0))/3, (9) 



where the pak operator defined as 

N 

PAK(t) = $:^^^^''-''^^*^=/^-AKW (10) 

8 



is the Fourier component of the adsorbate number-density operator, 

N 

p(R,t) = ^5(R-R,(t)). (11) 

i=i 

Thus, in the brackets denote the ensemble average over the trajectories associated with 
each adsorbate Rj(t). The intermediate scattering function is the typical observable issued 
from HeSE and NSE experimental techniques. Taking into account the relations ([7]) to 
(fTTj) . we note that the dynamic structure factor can be expressed in terms of a density- 
density correlation function and determined by the spectrum of the spontaneous fluctuations. 
Moreover, the static structure factor, defined as S'(AK,t = 0), is related to g(R), which 
describes the instantaneous correlation between adsorbates. 

From relations ([3]) and one finds by direct differentiation in the reciprocal space the 
continuity or number conservation equation, 

^^ = zAK.W), (12) 

which inserted into and taking into account (jTj), we obtain the basic relation 

SiAK,u) = ^J;(AK,a;), (13) 

where Ji is the longitudinal (projected along the wave vector transfer direction) current 
correlation function. This relation is rigorous as far as the adparticles are not being created 
or absorbed. 

The dynamic structure factor S'(AK, u) is an even function in frequency. Therefore, 
all the odd moments vanish and the frequency moments or frequency sum rules (static 
correlation functions) are AK-dependent quantities, 

'(92"/(AK,t)" 



w2"(AK) = — / w2"5(AK,a;)rf, 
2tt J 



u = -1 



in 



J t=0 



(14) 



which are obtained from a Taylor series expansion of the intermediate scattering function 
at short times (or small distances). In particular, the zeroth-order moment corresponds to 
the static structure factor, 

a;°(AK) = ^(AK), (15) 
since it describes the average distribution of interparticle distances on the surface; the second 



order moment is related to the thermal speed vq = \f^lk^Tjrn, as 

a;2(AK) = AK^t;^. (16) 



Due to the quantum character of the different operators introduced above, several com- 
ments are worth stressing. First, pAK^t) and P^k(O) commute only at t = 0. Second, the 
system studied here is assumed to be stationary and, therefore, the origin of time is arbitrary 
for the correlation function associated with the density operators. Third, the complex char- 
acter of the corresponding correlation function is a signature of the quantum dynamics of 
the interacting system. Fourth, the G-function is also complex, but the dynamic structure 
factor is real and positive definite because it represents a cross- sect ion. More properties 
of the pak(^) ojjerator, the G-function and the dynamic structure factor can be found in 
ovesey's book [^]. And fifth, the so-called detailed balance condition can be expressed as 



36| 



S{AK, u) = e'''^'^S{-AK, -u), (17) 



which expresses that the probability that a He atom loses an energy hu is equal to e'*^^ 
times the probability that a He atom gains an energy hw. 

After van Hove if -Ro is the range of the G-function and Tq its relaxation time, h/ Rq and 
H/Tq determine the orders of magnitude of average momentum and energy transfers in the 
scattering process of the probe particles, which for light masses display the observable recoil 
effect. Thus, the time variation of G affects the total scattering and angular distributions 
only for a particle spending at least a time of order Tq over a correlation length Rq. Moreover, 
if the mean de Broglie wavelength, A^, is small compared to interadparticle distances or the 
range of adsorbate-adsorbate interaction, no quantum effect will manifest in the G-function, 
which deals with pairs of adparticles separated by distances of the order of Rq. Nevertheless, 
for small timescales (t ^ Tq or t ~ the dynamics entirely concentrates in a region of 

the order of or less than Ab, and quantum effects are noticeable. At these distances, the 
adparticles can be considered as a two-dimensional free gas. The imaginary part of the 
G-f unction is greater at small values of time. 

B. The Hamiltonian for the system and the thermal bath 

1. The one-bath model 

In order to go a step further into the dynamics, we need to specify a Hamiltonian as 
introduced in Eq. ([8]). In surface diffusion, the full system-|-bath Hamiltonian is usually 
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written [l2| as 



H 



AT 



i=l 
N 

+E 



i=l 



2mi 2 

,2 



■ X 



2mi 



m 

2 



(18) 

where {Px,Py) and {x,y) are the adparticle momenta and positions, and {Pxn^i) and {py-,yi) 
with z = 1, ■ ■ ■ , are the momenta and positions of the bath oscillators (phonons), with 
mass and frequency given by rrii and Ui, respectively; phonons with polarization along the 
z-direction are not considered. The Hamiltonian was originally considered by Magalinskii 
[37 ] and Caldeira and Leggett [38], who used it for weak and strong dissipation (a general 
discussion about the Hamiltonian (fT8|) can be found in Weiss' book 39|). In surface dif- 
fusion, V{x, y) is in general a periodic function describing the surface corrugation at zero 
temperature. 

The harmonic frequencies of the bath modes and the coupling coefficients are expressed 
in terms of spectral densities, defined as 



N 



(19) 



with i = x,y. This enables the passage to a continuum model. The associated friction 
functions are defined through the cosine Fourier transform of the spectral densities as 

2 



= / d,^ cosut, 

vrm Jq u 



(20) 



with i = x,y. For Ohmic friction, f]i{t) = 2rii5{t) where rji is a constant and 5{t) is the 
Dirac delta function. In this model, the noise is shown to be white when Ohmic friction 
is assumed. The paradigm of white noise is the Gaussian white noise. Dealing with large 
systems (the surface seen as a thermal bath) where the number of collisions between substrate 
and adsorbate is very high, the fundamental theorem of probability theory, namely the 
central limit theorem, assures that the fiuctuations of the bath will be Gaussian distributed. 
Diffusion can then be described by a Brownian-type motion involving a continuous Gaussian 
stochastic process. In virtue of the FD theorem, such fiuctuations can be related to the 
friction coming mainly from surface phonons: the phonon friction. Electronic friction due 
to low-lying electron-hole pair excitations are usually neglected in most of cases. 
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2. The. two-bath model 



The diffuse elastic intensity of the He atoms scattered at large angles away from the 
specular direction provides very detailed information on the mobility of adsorbates on sur- 
faces. Based on the transition matrix formalism, Manson and Celli 10| proposed a quantum 



diffuse inelastic theory for small and intermediate coverages of adsorbates on the surface 
by ignoring multiple scattering effects of He atoms. The dynamical structure factor is then 
obtained by assuming all the crystal vibrational modes (A^) and point-like scattering centers 
(M) satisfying the harmonic approximation with a given frequency distribution function. 
Therefore, following the same type of arguments, we could assume two independent, un- 
correlated baths to describe diffusion of interacting adsorbates. As before, the first bath 
consists of the surface phonons. Meanwhile, the second bath is formed by M adsorbates 
which, obviously, changes with the surface coverage given by experimental conditions 28 1. 

For a two-bath model (for a given coverage), we take one adsorbate as the tagged particle 
or system, while the remaining ones constitute the second bath descibed by M harmonic 



oscillators. In this way, the corresponding total Hamiltonian reads as 

H 



28| 



2m 2m 



N 



i=l 

N 

+E 

i=l 

M 

M 



2m.; 2 



2m, 
2m 



mi 



+ 



X 



p: 

2m 



d 



y 



(21) 



where the barred magnitudes label the same quantities as in the one-bath model, but now 
referring to a bath of M adsorbates, which are also considered as a collection of harmonic 
oscillators. The c^^. and (i^,^, coefficients, with k = x,y, give the coupling strengths between 
the adsorbate and the substrate phonons or other adsorbates, respectively. The spectral 
density for the two baths is defined analogously to the one-bath model. 



Ji(oj) 



N 

-y 

2 ^ 



[5(u;-u;g]+-^ 



M 



h 



(22) 
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but now it is split into two sums, one spectral density due to the surface phonons and the 
other one due to adsorbates. In a similar way, the friction functions are defined as in Eq. fl20|) 
but now the spectral density is given by f l22p . The friction function is also split into two 
terms, one due to the phonons, 7(t), and the other one due to the presence of adsorbates, 
A(t). This last friction function could be interpreted like a collisional friction. Again, for 
Ohmic friction, T]i{t) = 2(7 + X)6{t), where 7 and A are constant and 6{t) is the Dirac delta 
function. 

Extension of the one-bath model to two baths is carried out to describe collisions among 
adsorbates when the coverage is increased up to a certain value. Again, if the friction 
is assumed to be Ohmic, the noise will be white. Adsorbate collisions can be seen as 
discrete events. It is well known that an appropriate way to model this type of noise is by 
considering a Poisson process, being generally designated as a white shot noise. As has been 



shown elsewhere 



21 



22l |. this white shot noise can be obtained as a limiting case of a color 



noise. It is also clear that when the diffusion regime is reached, the discrete Poisson process 
becomes a continuous Gaussian process. The introduction of the second bath allows us to 
describe adsorbate collisions by a new friction coefficient, the collisional friction. If the two 
baths are not correlated, the corresponding noises are also not correlated and this is the key 
point of the ISA model. The total friction coefficient, the sum of the phonon and collisional 
frictions, is related to the fluctuations of both baths through the FD theorem. This sum of 
frictions has recently been estimated by the He spin echo technique 27|, l31 |. 

Before concluding this subsection, we would like to emphasize that if an external driving 
force is also present (He atoms are sampling the surface to explore the motion of adparticles) 
a new extra term should be added to the total Hamiltonian, (ITSl) or (12T]) (see Section UTTl) . 



C. The Langevin equation in the two-bath model 

The next task is to eliminate the environmental (two-baths) degrees of freedom, which 
leads to a damped equation of motion of the system coordinates. If the Heisenberg picture 
of quantum mechanics is used, where the time evolution of a given operator A is given by 

A=Uh,A], (23) 
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a generalized Langevin equation (GLE) for each coordinate of the system is obtained, reading 
as follows 

mx{t) +m I 7],{t - t') x{t') dt' + ^}-^pyl = N^{t) (24) 



and 



my{t) +mj'^ r^y{t - t') y{t') dt' + ^^-^^ = Ny{t), (25) 

where the associated friction functions are defined through the cosine Fourier transform of 
the spectral densities given by Eq. fl22l) . 

2 /"^ , JAio) , , 

riAt) = / duj coswt, (26) 

vrm Jo w 



with i = x,y. Due to the splitting of the spectral density (|22|) . note that the friction function 
also splits into two terms, one due to the phonons, 7(t), and the other one due to the presence 
of adsorbates, \{t). 

The inhomogenity in Eqs. and ( 12 5 p represents a fluctuating force which depends 
on the initial position of the system and initial positions and momenta of the oscillators 
of each bath; a generalization of the one bath model 39|]. The fluctuating force in each 



direction can again be split into two sums, one due to the phonons and the other due to 
adsorbates. As both baths are assumed to be uncorrelated, the same property holds for the 
two noises. For each noise and each cartesian component, it can be easily shown that its 
equilibrium (canonical ensemble) expectation value with respect to the heat bath including 
the corresponding bilinear coupling to the system vanishes. On the contrary, the noise 
autocorrelation function (each cartesian component and each bath) is a complex quantity 
because in general it does not commute at different times. In the classical limit /?. — ?■ 0, 
each noise correlation reduces to mkBTr]i{t), with i = x,y. For Ohmic friction, i.e., delta 



correlated, we have white noises. In the quantum case 39 



40[, and also for Ohmic friction. 



the imaginary part of each noise function is a step function and its real part goes with 
csch^(7rt/^/3). Thus, at zero surface temperature, the noise is still correlated even for long 
time (it decays algebraically like t~^) in contrast to the classical case. These facts give rise 
to important differences with respect to the classical case such as, for example, the noise and 
the system coordinates are correlated instead of being zero. In order to simplify the theory, 
we will consider only classical noise but keeping in mind that our quantum results (see 
Section IIIip will be only valid for not too low surface temperatures. In the QHAS and HeSE 
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experimental techniques used for fast diffusion, the lowest attainable surface temperature 
is around 50-100 K. In a future work, quantum noise and tunneling will be considered in 
surface diffusion problems at very low temperatures. 

Thus, if 'r]{t) = 2ri5{t) = 2(7 + X)6{t) (Ohmic friction), Eqs. and ([25]) reduce to two 
coupled standard Langevin equations (Markovian approximation) (the delta function counts 
only one a half when the integration is carried out from zero to infinity) 



20I-22 



221, 

mR= -m?7R-F(R),+5N, (27) 

which is the basis of the ISA model, and where (5N is given by the sum of two noncorrelated 
noises: the lattice (thermal) vibrational effects and the adsorbate-adsorbate collisions, which 
are simulated by a Gaussian white noise (G) and a shot white noise (S) (Poissonian dis- 
tributed), respectively. Thus, for each degree of freedom, we have 6N{t) = dNcit) + 6Ns{t). 
Let us remark that a Poissonian distribution behaves as a Gaussian distribution for very 
long times. The Langevin equation is then solved for one single particle in presence of two 
noises. 

For a good simulation of a diffusion process, one has to consider very long times in 
comparison to the timescales associated with the friction caused by the surface or to the 
typical vibrational frequencies observed when the adsorbates keep moving inside a surface 
well. This means that there will be a considerably large number of collisions during the 
time elapsed in the propagation, and therefore that, at some point, the past history of the 
adsorbate could be irrelevant regarding the properties we are interested in. This memory 
loss is a signature of a Markovian dynamical regime, where adsorbates have reached what 
we call the statistical limit. Otherwise, for timescales relatively short, the interaction is 
not Markovian and it is very important to take into account the effects of the interactions 
on the particle and its dynamics (memory effects). The diffusion of a single adsorbate 
is thus modeled by a series of random pulses within a Markovian regime (i.e., pulses of 
relatively short duration in comparison with the system relaxation) simulating collisions 
among adsorbates. In particular, we describe these adsorbate-adsorbate collisions by means 
of a white shot noise as a limiting case of a colored shot noise |25|, as mentioned above. This 
interaction is, therefore, described in terms of the collisional friction, which depends on the 
surface coverage. Thus, the ISA model essentially consists of solving the standard Langevin 
equation with two noise sources and frictions: a Gaussian white noise accounting for the 
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friction with the substrate and a white shot noise characterized by a colhsional friction 
simulating the adsorbate-adsorbate colhsions. This way of simulating the interaction among 
adsorbates reduces the dynamical problem to the diffusion of a single adsorbate (like in the 
SA approximation) and the factor appearing in the S or I functions, Eqs. ([7]) and (Q, 
has not to be considered. 

Finally, the surface coverage 6 and A or collisional friction can be related in a simple 
manner. In the elementary kinetic theory of transport in gases ^ diffusion is proportional 
to the mean free path /, which is proportionally inverse to both the density of gas particles 
and the effective area of collision when a hard-sphere model is assumed. For two-dimensional 
collisions the effective area is replaced by an effective length (twice the radius p of the 
adparticle) and the gas density by the surface density a. Accordingly, the mean free path is 
given by 

r=— (28) 

Taking into account the Chapman- Enskog theory for hard spheres, the self-diffusion coeffi- 
cient can be written as 

D = — \[^. (29) 

Qpa V m 

Now, from Einstein's relation, and taking into account that 9 = a^a for a square surface 
lattice of unit cell length a, we finally obtain 



A = 5fJ^«. (30) 

Therefore, given a certain surface coverage and temperature, A can be readily estimated from 
fl30|) . Notice that when the coverage is increased by one order of magnitude, the same holds 
for A at a given temperature. Notice that A could also be considered as a phenomenological 
parameter, just like the substrate friction 7. 

D. Linear response functions 

The dynamic structure factor can also be related to the linear response function of the 



system 



6| 



0(AK,t) = -^([pAK(t),pkK]), (31) 
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through the FD theorem, expressed as 

S{AK, uj) = [1 + n{u)] / e"^' 0( AK, t) dt, (32) 

where 1 + n{uj) = [1 — exp(— /iw/?)]"^ is the detailed balance factor, with n{uj) the Boltzman 
factor. Equation ( l32ll relates the spectrum of spontaneous fluctuations, ^(AK,^;), to the 
dissipation part of the response function. Time derivatives of the response function are 
related to the Heisenberg equation of motion fl23|) . of the pak operator and moments of the 
dynamic structure factor involve nested commutators to evaluate them. The 0-function is 
a causal function because it can not be defined before the external perturbation has been 
switched on. For the scattering with He atoms, the perturbation starts at — oo and ends at 
+00 and typically has a bell shape (see Section UTTl) . 

In ( l32l) . the time Fourier transform of defines a generalized susceptibility function 
x(AK,co') and, therefore, can be reexpressed as 

S{AK, uj) = -i[l + n(w)]x(AK, uo). (33) 

This susceptibility is complex and the real and imaginary parts are related through the well- 
known Kramers-Kronig or dispersion relations [^. On the other hand, the time derivative 
of the linear response function is related to the so-called relaxation function as follows 

/oo 
c/tXAK,t'), (34) 

which describes the relaxation of the density after the external perturbation (He atoms) 
has been switched off. Physically, i?,^ goes to zero as t — oo and at t = this function 
coincides with the isothermal susceptibility. Thus, the dynamic structure factor can again 
be reexpressed in terms of the relaxation function as 

S{AK,uj) = [l + n{io)]ujR^{AK,uj), (35) 

where i?0(AK,a;) is the time Fourier transform of the relaxation function. 

Finally, the dynamic structure factor can also be expressed in terms of the Green function 
which is directly related to the linear response function, generalized susceptibility and the 



relaxation function 



6|, 



S{AK,u) = -\ [1 + n(w)]ImG(AK,a;), (36) 
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where 

ImG(AK,a;) = -hlmx{AK,uj), (37) 
with ImG'(AK,a;) being the imaginary part of the Green function. 

III. APPLICATIONS 

In general, an exact, direct calculation of /(AK,t) or ^(AK,^;) is difficult to carry out 
due to the noncommutativity of the adparticle position operators at different times, which 
obey a Langevin-Markovian equation, here described by ( 12 7p . where the friction is assumed 
Ohmic and the interaction potential is not separable. However, for certain simple cases, 
close formulas can be easily obtained. 

The product of the two exponential operators in ([9]) can be evaluated according to a 
special case of the Baker-Hausdorff theorem (the disentangling theorem). If A and B are two 
operators then e^e^ = ^{^M/'^^a+b ^ which only holds when the corresponding commutator 



is a c-number. Thus, (El) reads now as 



34| 



/(AK, t) = /i(AK, t)/2(AK, t), (38) 

which is a product of two quantum intermediate scattering functions Jj(AK, t) with j = 1,2 
associated with the exponentials of the commutator [A, B] and A + B, respectively. If we 
identify 3J] the operators A and B a.s A = iAK. ■ R(0) and B = — iAK ■ R(t), the factor 
Ji involving their commutator will depend on the character of the dynamics; for classical 
dynamics, this factor is one. The second factor can also be written as follows 

within the so-called Gaussian approximation and where C^,(^) = (f AK(i)'yAK(0)) is the veloc- 
ity autocorrelation function along the direction given by AK or the longitudinal direction. 
Equation ( l39l) is exact if the velocity operator gives rise to a Gaussian stochastic process. 
In the cases we are going to analyze below, we will discuss the factorization given by ( 138|) 
in more detail. 



18 



A. Diffusion on flat surfaces 



In the case of diffusion on flat or very low corrugated surfaces, where the role of the 
adiabatic adsorbate-substrate interaction potential is negligible and only the action of the 
thermal phonons and surrounding adsorbates are relevant, one can assume V{x,y) ~ 0. 
Thus, the stochastic single-particle trajectories R(t) running on the surface obey the follow- 
ing Langevin-Markovian equation ([2] 



mR(t) = -mT]'R{t) + 5N(t), (40) 

where (5N(t) = (5NG(t) + 5J>is{t) is the two-dimensional fluctuation of the total noise acting 
on the adparticle. 

He atoms are usually the probe particles and it is generally assumed that they do not 
influence the surface dynamics, that is, their influence can be considered a perturbation. 
Some considerations are in order about the driving force or external perturbation. First, the 
attractive part of the interaction potential does not play an important role in vibrational 
excitation problems, and it can usually be neglected. And second, the incident energy for 
the incoming particles is large compared to the vibrational excitation of the adsorbate (the 
frustrated translational or rotational modes which are of very low frequency). A semiclassical 
description of the projectile-adsorbate interaction allows for an estimate of the coUisional 
time and the duration of energy transfer. Thus, if r is the distance between the He atom 
and the center of mass of the adsorbate and their interaction is accepted to be exponentially 



repulsive 4l|, K(r) = y4exp[— a'r], it can be easily shown that the external force can be 
expressed as Fe{t) = Csech^a't with C = a'vi sin 9i/2 where Vi and 9i are the incident velocity 
and angle, respectively. The parameter a' gives the rate of energy exchange between the 
translational (He atoms) and frustrated translational (adsorbates) motions. The hyperbolic 
function has the physically correct behavior at the asymptotic limits {t — ?■ ±oo), and it 
is maximum at t = where the closest distance to the adsorbate is reached (bell shape). 
A similar expression can be obtained for the external force if instead of a pure repulsive 
interaction a Morse interaction is used. 

Due to the fact that for a flat surface no direction is priveleged, and if the adsorbate 
motion is driven by the external force Ff.[t), we have from Eq. (HU]) that 

{x{t))p + r^{x{t))p = ^F,{t). (41) 

lit 
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Within the hnear response framework we could write the particular solution of the differential 
equation as 

{S:{t))p= I ds(t){t-s)F,{s), (42) 

J — oo 

or, after Fourier transforming, 

{R{uj))p = xico)F,{co), (43) 

where 

F,(a;) = — CSC — . (44) 



a'2 \2a' 

Now, if we assume that cj/a' keeps close to unity during the interaction along time, the 
cosecant function will be also close to unity and the dynamic susceptibility will be that of a 
free adparticle on a flat surface, 

X(^) = - / . , (45) 

which is exact whenever an Ohmic friction r] is assumed and any direction given by AK is 
considered. This expression of the dynamic susceptibility is valid for both the classical and 
quantum case. 



1. Classical dynamics 



The adparticle motion can then be regarded as quasi-free since it is not ruled by a 
potential, but only influenced by the two stochastic forces. Within this regime, the velocity 
is a Gaussian stochastic process and the velocity autocorrelation function in any direction 
(since there is no priveleged direction) is given by Doob's theorem 22| . 



-r]t 



m 



(46) 



The expression for the intermediate scattering function resulting from (139!) . where 
I{AK,t) = hi^K^t) is 



/(Air, t) = exp [-x^ (e-^* + r]t - l)] 



(47) 



where the so-called shape parameter x 2a, |42| is defined as 



(48) 
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From this relation we can obtain both the mean free path I = t^J {vq) and the self-diffusion 
coefficient D = t{vq) [Einstein's relation). When the coverage increases, the collisions 
among adsorbates are also expected to increase, and so A and therefore rj. As can be easily 
shown, fj47j) displays a Gaussian decay at short times that does not depend on the particular 
value of 7] (ballistic motion), while at longer times it decays exponentially with a rate given 
by ri~^ . Thus, with r/, the decay of the intermediate scattering function becomes slower. 

The above described effects can be quantitatively understood by means of the expression 
of the dynamic structure factor obtained analytically from P7|) . 

711] I ) 

TT ^ n! u'^ + [{x^ + n)riY ' 

Here, the F symbol in the first line denotes both the Gamma and incomplete Gamma 
functions (depending on the corresponding argument), respectively. As can be noted in 
the high friction limit, fH9|) becomes a Lorentzian function, its full width at half maximum 
(FWHM) being F = 2?7x^, which approaches zero as rj increases (narrowing effect). This is 
in sharp contrast to what one could expect — as the frequency between successive collisions 
increases one would expect that the line shape gets broader (effect of the pressure in the 
spectral lines of gases). The physical reason for this effect could be explained as follows. 
As 7] increases particle's mean free path decrease and therefore correlations are lost more 
slowly. In the limit case where friction is such that the particle remains in the same place, 
the van Hove function becomes a 5-function, the intermediate scattering function remains 
equal to one and the dynamic structure factor consists of a 5-function at u; = 0. Conversely, 
in the low friction limit the line shape is given by a Gaussian function, whose width is 



F = 2\/2 \n2^ykBT/m AK, which does not depend on rj. This is the case for a two- 



dimensional free gas 



43 



44| . This gradual change of the line shapes as a function of the 



friction and/ or the parallel momentum transfer leading to a change of the shape parameter 

nn 

X is known as the motional narrowing effect [15|, 116|. Notice that, in our approach, friction 
is related to the coverage. Thus, at higher coverages a narrowing effect is predicted for a 
flat surface [2C 
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2. Quantum dynamics 



In the Heisenberg representation, Eq. (HOl) still holds, its formal solution being 

R(t) = R(0) + <^(rit) + / ^(rit - rit')SN(t')dt' , (50) 

mr] mr] Jq 

where P(0) is the initial adparticle momentum operator and $(a;) = 1 — . Then, the 
commutator between R(0) and R(t) is obtained from Eq. (!50l) obtaining a c-number. Then, 
the factor Ji, by assuming a classical noise as previously mentioned, reads as 34 1 



/i(AK,t) = exp 



<l>(r/t) 



exp 



h 7] 



(51) 



2rjm 

where Er = /i^AK^/2m is the adsorbate recoil energy. As is apparent, the argument of the 
exponential function becomes less important as the adparticle mass and the total friction 
increase. The time-dependence only comes from ^{rjt). At short times (< h(3), ^{r]t) ~ r]t 
and the argument of Ji becomes independent of the total friction, thus increasing linearly 
with time. On the other hand, in the asymptotic time limit, this argument approaches a 
constant phase. 

In order to calculate the I2 factor we start from Eq. (I40l) describing an adparticle with 
mass m moving on a flat surface in presence of more adsorbates. The dynamic susceptibility 
is also given by Eq. fH5|) and its time behavior by 

X(t) = —e-'''/^smh(r]t/2)Q(t) (52) 
mrj 

where B(t) is the step function due to causality. The FD theorem allows us to express the 
equilibrium position autocorrelation function, Cx{t) = (x(t)x(O)), in terms of the imaginary 
part of the dynamic susceptibility and, after Fourier transforming, we find 

CM) = / duo , ' ^ 53) 

From the relations 
and 

o / °° 2 \ 

coth(^M2) = ^ 1 + 25:-^ . (55) 



where the so-called Matsubara frequencies are defined by 



= ^, (56) 
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we can decompose the correlation function as Cx(t) = Sx(t) +iAx(t), i.e., into its symmetric 
and antisymmetric parts, respectively. For t > 0, these functions read as 

SAt) = ftsign[t] + ^e-A +^f2 (57) 

and 

Ax{t) = -^{l-e'^'), (58) 

which can be trivially related to ( l52l) by the FD theorem. In ( 1571) . the sign function of the 
real number t is defined as follows: +1 for t > and — 1 for t < 0. Now, since 

C.{t) = -^Cx{t), (59) 

then 

C.W ^ - a) e-.--|Lf: "r"' . (60) 

V am 2m I pm ^-^ m — vi, 

\' / n=l " 

where the real part is identical to H6|) except for the infinite sum of the Matsubara frequencies. 
Quantum effects are important at low surface temperatures, the long time behavior being 
mainly determined by the first term of the Matsubara series. In such cases, relaxation is no 



longer governed only by the damping constant 



39 



4Q| 



Finally, (l60ll is substituted into ( l39l) in order to obtain the factor I2, 

J2(AK,t) = Q-^K^ !Q{i-i')Cv{t')dt' ^ g-Aif2[/(t)+g(t)]^ ^g^^ 



where the time-dependent functions f{t) and g{t) are given by 

mf3ri'^ 2mrj J 



fit) = ( ^ - ^7^] [e'"* + Vt-l] (62) 



and 



The total intermediate scattering function (138!) can be then expressed as 

/(AK, t) = e-x'[-*^*-*(';*)]e-^^'^'W, (64) 

with X = ^K'^i'^o) /v"^ a = 1 + ihf3/2, the thermal square velocity being {vq) = l/m/3. 
The recoil energy is included in the imaginary part of the product x^Q^*? which disappears 
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when h ^ 0. Equation (16T]) is slightly different to that obtained elsewhere [3J], where the 
factor I2 was treated classically, 

/2(AK,t) =e-'^'[''*-*(''*)l, (65) 

and, therefore, the total intermediate scattering function fl38l) can be expressed as 

/(AK, t) = e"^' e"^' . (66) 

Equation fl64l) is the generalization of the intermediate scattering function for the quantum 
motion of interacting adsorbates in a fiat surface. The dependence of this function on 



AK^ through the shape parameter is the same as in the classical theory 22j. No previous 
information about the velocity autocorrelation function is needed. However, classically, the 
intermediate scattering function is usually obtained from Doob's theorem, which states that 
the velocity autocorrelation function for a Gaussian, Markovian stationary process decays 



exponentially with time [45|]. The ballistic or free-diffusion regime and the diffusive regime 
are apparent from fl66|) . The first one is dominant at very low times, rjt <^ 1, and the second 
one at very long times, rjt^l. 

The diffusion coefficient can be obtained from 

D = lim j C^{t')dt'. (67) 







Thus, from Eqs. fl60l) and fl62|) . the diffusion coefficient is the real part of the complex number 
given by 

D = ^-i—, 68 
mrj 2m 

which coincides with Einstein's law for the classical case (insuring that the adparticle velocity 
distribution becomes Maxwellian asymptotically). The same result is reached from the mean 
square displacement, ((R(t) — R(0))^), which takes into account only the symmetric part 
of the position autocorrelation function. Quantum fluctuations (the Matsubara frequencies) 
do not affect this result at low temperatures except the time limit to which the mean square 
displacement (MSD) is linear with time may become very large. At zero temperature, 
D is also zero and the MSD is no longer linear with time. The infinite sum of Matsubara 
frequencies determines now the long time limit behavior. As previously mentioned, the limit 
to zero surface temperatures is questionable if in the commutator we neglect the correlation 
between the noise and the coordinate system. 
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B. Diffusion on harmonic potentials 



The harmonic model is an appropriate working model to understand the bound motion 
inside the wells of a corrugated surface and, therefore, to also understand the behavior 
associated to the T-mode, which comes precisely from the oscillating behavior undergone 
by the adparticle when the diffusional motion is temporarily frustrated. Now if we again 
assume that lojoi keeps close to unit during the interaction along time, the cosecant function 
will be close to one and the dynamic susceptibility will be that of an adparticle subject to 
a one-dimensional harmonic potential 

= - 2 ^ ^ 2 > (69) 

Uq being the frequency of the harmonic oscillator. This expression is exact whenever an 
Ohmic friction rj is assumed. This expression of the dynamic susceptibility is valid for both 
the classical and quantum cases. 



1. Classical dynamics 



In contrast with the case of a dynamics where V{x, y) does not play a relevant role, 
we can devise a particle fully trapped within a harmonic potential well. Thus, for a har- 
monic oscillator, the behavior of the adparticle becomes very apparent when looking at the 



corresponding velocity autocorrelation function, which reads 



15 



22,|45| as 



with 



m uo 



UJ 



-r,t/2 



COS (cut + 5i 



T 



(70) 



(71) 



and tan Si = r]/2u). Note that ( H6l) can be easily recovered after some algebra in the limit 
Uo^O from (1701). 

The only information about the structure of the lattice is found in the shape parameter 
through AK [see Eq. fHHl) ]. When large parallel momentum transfers are considered, both 
the periodicity and the structure of the surface have to be taken into account. Consequently, 
the shape parameter should be changed for different lattices. The simplest model including 
the periodicity of the surface is that developed by Chudley and Elliott [461], who proposed a 
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master equation for the pair-distribution function in space and time assuming instantaneous 
discrete jumps on a two-dimensional Bravais lattice. Very recently, a generalized shape 



parameter based on that model has been proposed to be [42 1 



where, within our approach, 7 is substituted by rj. Here, ry(AK) represents the inverse of 
the correlation time and is expressed as 

r,(AK) = zy ^ Pj [1 - cos(j ■ AK)], (73) 
j 

V being the total jump rate out of an adsorption site and Pj the relative probability that a 
jump with a displacement vector j occurs. 

Substituting now flTOj) into fl39|) leads to the following expression for the intermediate 
scattering function 

I{AK, t) = exp <^ [cos 5i - 6-"*/^ cos{cut -5i)]\. (74) 

The argument of this function displays an oscillatory behavior around a certain value with 
the amplitude of the oscillations being exponentially damped. This translates into an also 
decreasing behavior of the intermediate scattering function, which also displays oscillations 
around the asymptotic value. This means that after relaxation the intermediate scattering 
function has not fully decayed to zero unlike the case of absence of a potential. Again, in 
the limit Uq ^ 0, (!7i|) approaches (1471) . 

In order to obtain an analytical expression for the dynamic structure factor, it is conve- 
nient to express fl7i|) as [22] 

°^ / 1 \ m+n 

= g-X?^i ^ ^2(m+n)^m^ngi(m-n)5ig-(m+n)?7V2gj(m-n)(I.t^ ^j^'j 

m,n=0 

where 
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with 

^ UJQ 7]^{2{r)/2)OsmSi + [u^ - (77/2)^] cos(5i} 
' CO [{r]/2y + u^f ' ^ 

^3 = / (78) 

" io2{r]/2 + tuy' ^^^^ 
where the coefficients Ai has been put in terms of rj, u and 61. From (1751) . it is now 
straightforward to derive an expression for the dynamic scattering factor, which results 

vr m! n! 

r?i,n=0 

XT 1 ^^V^M^ rvi^. (80) 



For a harmonic oscillator, there is no diffusion and, therefore, flHOj) is only valid when 
m ^ n. All the Lorentzian functions contributing to ( IHOl) are due to the creation and 
annihilation events of the T mode. These Lorentzians are characterized by a width given by 
r = (m + n)ri/2, which increases with 7]. This broadening (proportional to r/) undergone by 
the dynamic structure factor is thus contrary to the narrowing effect observed in the case of 
a fiat surface. It can be assigned to the confined or bound motion displayed by the particle 
ensemble when trapped inside the potential wells. Hence, in order to detect broadening of 
the line shapes in surface diffusion experiments, adparticles must spend some time confined 
inside potential wells, since the broadening will be induced by the presence of temporary 
vibrational motions. 



2. Quantum dynamics 

The formal solution of Eq. (1271) is given by 

R(t) = R(0) + ^ <!>{r]t) + — f - vt')[F{W)) + SN{t')]dt', (81) 
mrj mrj Jq 

where the force F is given by Hooke's law and P(0) is the initial adparticle momentum 
operator and $(x) = 1 — e~^. The presence of the adiabatic force introduces an additional 
commutator, [Ro,F(R(t))] = (z/i)(9F(R(t))/9Po, where the dependence of the adiabatic 
force on the initial state (R(0), P(0)) is through R(t) which is again negligible in a quantum 
Markovian framework [3J]. Thus the factor Ji is the same as for a fiat surface, Eq. (|5T|) . 
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In order to calculate the I2 factor we again start from Eq. (1271) . The dynamic susceptibility 
is also given by f l69|) and its time behavior by 

Xit) = —e-"'/^ smcoteU), (82) 

where Q(t) is again the step function due to causality. According to the FD theorem, as 
before, the equilibrium position autocorrelation function can be expressed as 

CM) = / du— ^ — ^ (83) 

whose symmetric and antisymmetric parts are given, for t > 0, as 

Sxit) = — 7^ n P COS tut + h/2) sin cot] - —— > — (84) 

m[5uul^ ' ^ /3m ^ (77/2 - z/„)2 + £i;2 ^ ' 

and 

A, it) = —e-'''/^ sin ut, (85) 

2mijj 

respectively, and the velocity autocorrelation function by 

Q 00 2 —u t 'fe 2 

C.{t) = -^e-^^/'cos{ut + S,) + ^Y. f /o"^"' ^2 ^, -'f^e-^^/'cos{ut + S,), (86) 
mpu pm [rj/Z — Unj — w 2mu 

with tan52 = rju/ujQ. Again, the real part is the same as in the classical case except the 
Matsubara series. The same considerations about the surface temperature in the quantum 
regime can be done as before. The I2 factor can be expressed as in (|6T1) through the functions 

i 9 ,,2 

m Ml,! nm < ^ 



mPuj f3m (77/2 — — uj"^ 

iht 



2m 



2uJqUJ 



(87) 



and 



g(t) = {d; -e~^*/^[a;cosa;t + (a;n^ + 7/2)sina;t]| 



2 ^ 1 - e ''"*(z/„t + 1) ihwn , , 



/3m (77/2 — z/„)2 — ci;2 2muj 



with 



^0 = 7-t(^' + 2^^o^' - 2a;o'7 - 4a;o), (89) 



(7i(t) = — ^[uUq{uq — T] )t + 2ujQriuj — rj u] cos ut, (90) 





-Vt/2 



92it) = ^^[ir]/2)oj',i3iu', - r]')t + Auj^ - uj^ - r/72] sincDt. (91) 
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The total intermediate scattering function will be the product of the factors Ji and I2 given 
by fl5T]) and fl6T|) . taking into account fl87|) and fl88l) . respectively. 



C. Diffusion on corrugated surfaces 



-Z. Classical dynamics 



The broadening of the diffusion line shapes has been shown to be produced by the tempo- 
rary trapping 20|. We start this subsection by considering a general velocity autocorrelation 
function which keeps the same functional form of Eq. ( I70|) . but whose parameters do not 



hold the same relations as those characterizing a harmonic oscillator 



15|, that is, 



m 



e"''* cos{ujt + 6), 



(92) 



where the values of the parameters fj, u and 6 are obtained from a fitting to the numerical re- 
sults issued from solving the standard Langevin equation with periodic boundary conditions 

one easily reaches the corresponding expression for the intermediate 



(EZD. From Eq. (|92 
scattering function, 



22| 



I{AK,t) = e-^'^'('^'*) 



E 



m,n=0 

which is analogous to fl75p . and where 



^2{m+n) ^mj^n^i{m—n)5 ^—{m+n)r}t^i{rn—n)uit (93) 



f{u, t) = Ai + A2t + A^/^e-^^-'^^* + Aie-''^e-^^+' 



and 



Ai 
A, 



rf [2rjuj sin 5 + (w^ — if ) cos 5) 



r,2^2 



{rf + Cj 
Tfijj CO?, 5 — Cjsm.5) 



(94) 

(95) 
(96) 
(97) 

2{fi + iQf ^^^^ 
Unlike the case of the harmonic oscillator, notice now that there is a new linear dependence 
on time in / because of the parameters fj, u and 5 are time-independent. This leads to an 



rf + a;^ 



rf 



2(r] — iuy ' 



rf 
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exponentially decaying factor in (p3|) . which accounts for the diffusion and that makes the 
intermediate scattering function to vanish at asymptotic times. In this sense, the intermedi- 
ate scattering function can be considered as containing both phenomena, diffusion and low 
vibrational motions. This effect is better appreciated in the dynamic structure factor, 

TT ^-^ mini 

Tn,n=0 

^ X?^2 + (m + n)fj ^gg^ 

[cu — (m — n)u]'^ + [x^^2 + {m + n)f/]2 

This general expression clearly shows that both motions (diffusion and oscillatory) cannot 

be separated at all. The Q-peak is formed by contributions where m = n, for which each 

partial FWHM is given by 

TQ = x^A2 + 2mf]. (100) 

Analogously, the T-mode peaks come from the sums with n ^ m and each partial FWHM 
is given by 

TT = xf^2 + {m + n)f]. (101) 

If the Gaussian approximation is good enough, the value of f] will not be too different from 
the nominal value of rj and, therefore, both peaks will display broadening as 1] increases. 
This is a very remarkable result since a relatively simple model, as the one described here, 
can explain the corresponding experimental broadenings observed with coverage. Thus, 
broadening arises from the temporary confinement of the adparticles inside potential wells 
along their motion on the surface 20|]. The problem of the experimental deconvolution has 
been already discussed elsewhere 3]. Here we would like only to mention that using this 
simple model, such deconvolutions would be more appropriate in order to extract useful 
information about diffusion constants and jump mechanisms. Finally, as mentioned before, 
the motional narrowing effect will govern the gradual change of the whole line shape as a 
function of the friction or, equivalently, the coverage, the parallel momentum transfer and 
the jump mechanism. 



2. Classical dynamics. The Kramer's turnover framework 



The theory of activated surface diffusion in one dimension was developed 47|, |48|] from 
Kramers' solution to the problem of escape from a metastable well. 49|, |50|] The underlying 
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dynamics is assumed well described by the Langevin equation provided that the reduced 
barrier height is of the order of ~ 3 or higher, the energy loss to the bath of trajectories 
close to the barrier top is given by classical mechanics and the potential at the barrier top is 
approximately parabolic. It has been shown that Kramers' based theory with finite barrier 
correction terms can then be replaced by Langevin numerical simulations 



The starting point is a kinetic equation for the stationary flux of particles exiting each 
well at either barrier. This flux is affected by the rate of particles exiting the jth well and 
those arriving at the well from the two neighboring wells, j + 1 and j — 1. Here we are going 



to give the main ana^ 
found elsewhere. 18 



ytical expressions derived from Kramers' theory, more details can be 



48| A central quantity in the theory is the reduced average energy 
loss 6 to the bath as the adatom traverses from one barrier to the next. For a single cosine 
potential as that considered here with barrier height = 2Vo, the energy loss is given by 

S^^, (102, 

where ooq = 27r^yVo/ma'^ is the harmonic frequency of oscillation near the well bottom, m 
is the mass of the adatoms and a is the unit cell length. Since typically many experiments 
or calculations are carried out under conditions of large reduced barrier heights, 6 can be 
unity or even larger, even though the damping constant is rather small. 

In the moderate to strong friction limit where the rate limiting step is spatial diffusion 
(sd) across the barrier, the rate of the escape from the well in both directions is given by 



the Kramers-Grote-Hynes formula 49|, |51| 



T,,= -^^exp[-VykBT], (103) 



where the Kramers-Grote-Hynes prefactor is 



which has been generalized to two baths and appears as a renormalization taking into account 
recrossings, since we are working implicitly in normal mode coordinates for the diffusing 
particle and the two baths. Finally, for the partial rates one finds 

rp27T I try PIT Jl 



r, = -— / rfAir sin^ (Air/2) cos{jAK) exp i - / dx In 

Jo \^ Jo 



1 - P\x) 
1 + P2(a;) cos(Air) 



(105) 
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where the function P{x) is defined as 

P{x) = exp 



6 



and AK is written dimensionless. The rate of escape from the zeroth well, 



(106) 



(107) 



and the relative probability or a jump of length j is given by the probability of being trapped 
at the jth well, 

P, = ^. (108) 

K 

For a one-dimensional periodic potential, the diffusion coefficient is related to the escape 
rate by 

1 1 °° 

Z} = -«:(a2) = -a2 5^j%, (109) 



where (a^) is the mean square path length. The diffusion coefficient can then be expressed 
in close form as 



D = DsdT-^ exp <( - 



TT 



7r/2 



dx\n[l + P{x)] 



(110) 



with = (l/2a^)rsrf the diffusion coefficient in the spatial diffusion regime, and T the 
depopulation factor for the metastable well first given by Melnikov j47| . 



T = exp 



TC 



J dx\n[l-P{x)]\ . 



Ill) 



In analogy to the Chudley-Elliott model 
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46| , an analytical expression for the FWHM 



of the dynamic structure factor can also be obtained by imposing a master equation for the 
intermediate scattering function. One easily sees that if the dynamic structure factor has 
the ubiquitous Lorentzian shape, the FWHM is given by 



Tr/2 



dx In 



P\x) 



(112) 



1 + P2(a;) _ 2P{x) cos(AK)_ 

This equation is important in the sense that assuming the validity of Kramers' model and 
the master equation approach, it allows for a direct comparison with the experimental data 
and/or Langevin numerical simulations and therefore an estimation of the spatial diffusion 
rate F^^ and the energy loss 5. From these parameters and their temperature dependence, 
one can further infer the barrier height, the friction coefficient and the barrier frequency. 
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In order to solve Eq. (127|) we have used the velo city- Ver let algorithm, which is commonly 



applied when dealing with stochastic differential equations [52|]. For the average calculations 
shown here a number of 10,000 trajectories is sufficient for convergence along the [100] 
direction. The initial conditions are chosen such that the velocities are distributed according 
to a Maxwell-Boltzmann velocity distribution at a temperature T (i.e., the ensemble is 
initially thermalized) , and the positions such that they cover the extension of a single unit 
cell of the potential model used (see below). Regarding the dynamical parameters, we have 
used the same parameters previously chosen for the Na atom on a Cu(OOl) surface, that 
is, 7 = O.Iuq, where Uq is obtained from the harmonic approximation assumed near the 
well bottom with a barrier of Vq = 41.4 meV. The mass and radius of the adparticles are 
those of a Na atom, since this adsorbate has been widely used in QHAS experiments over 
the Cu(OOl) surface (non-separable potential) [23|. As for the coverage, 9=1 corresponds 
to one Na atom per Cu(OOl) surface atom or, equivalently, a = 1.53 x 10^^ atom/cm^; 
a = 2.557 A is the unit cell length; and p = 2 A has been used for the atomic radius. Once 
the surface temperature and the coverage are fixed, the corresponding A value is obtained 
from Eq. ([30]). 

In Fig. [1] jump rates calculated in ps^^ at 200 K are shown as a function of the surface 
friction for a nonseparable interaction potential 23| and two coverages (solid line 6 = 0.028 
and dashed line, 6 = 0.18). These jump rates are obtained from Kramer's theory and squares 
and circles from mean first passage time calculations. As can clearly seen, the left part of 
the turnover region is slightly shifted depending on the coverage used. Notice that without 
including finite barrier corrections, the agreement is fairly good indicating that Kramers 
one dimensional theory can be convenient for interpreting QHAS measurements even when 
interacting adsorbates are present. 

In Fig. [2]the tracer diffusion coefficient in a.u. versus the total friction at 200 K is plotted. 
White circles are issued from the ISA model and black circles from Kramer's results. The 
agreement is again quite good and Einstein's law is clearly fulfilled. 

Finally, in Fig. [3] the FWHM (in /leV) of the Q-peak is shown for two different coverages 
(low, 0.028, and moderate, 0.18) at 200 K as a function of the parallel wave vector transfer 
covering the first Brilloiun zone. Results from the ISA model are plotted in black squares 
and circles. White squares and circles are experimental results and the solid lines are coming 
from Kramer's theory which are also obtained by assuming a Chudley-Elliott model. The 
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FIG. 1: Jump rates in ps~^ as a function of the coverage (solid line 6 = 0.028; and dashed line, 
6 = 0.18) for the diffusion of Na atoms on Cu(OOl) surface at 200 K, along the [100]. Kramers 

results are plotted in lines and mean first passage time calculations are plotted in circles and 
squares. 




2.4e-05 2.Se-05 3.2e-05 3.6e-05 4e-05 4.4e-05 4.8e-05 
Ti(a.u) 



FIG. 2: Tracer diffusion coefficient in a. u. as a function of the coverage for the diffusion of Na 
atoms on Cu(OOl) surface at 200 K, along the [100]. Kramers results are given by black circles and 
Langevin numerical simulations by white circles. 



agreement among theoretical results is fairly good but with the experimental results are 
poorer for the high coverage. Values of the coverage around 0.16 mark the upper limit 
where the ISA model can work. 
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FIG. 3: Full width at half maximum (in /ueV) of the quasielastic peak as a function of the parallel 
wave vector transfer for the diffusion of Na atoms on Cu(OOl) surface at 200 K and two coverages: 
0.028 (low coverage) and 0.18 (intermediate coverage), along the [100]. 

3. Quantum dynamics 

Our starting point is again Eq. ( 1271) where now the adiabatic force is derived by any 
general interaction potential. The same discussion as in the harmonic case can be followed 
for the Ji factor. 

For Na atoms, the pairwise interaction potential is repulsive and the mean interparticle 
distance should be most of the time greater than A^- Thus, the I2 factor could be replaced, 
in a first approximation, by the classical counterpart, Eq. ( |93l) . The error comes from 
small times but due to the fact the diffusion process is a long time one, the influence on the 
quasielastic peak (wave-vector dependence) and quantum diffusion constant (Einstein's law) 



will be really small for massive particles jSJ] . In Fig. HJ we show the effects of the quantum 
correction in the diffusion process studied here at two different surface temperatures for a 
coverage of 0.028. For comparison, in this plot, the classical intermediate scattering function 
and the real part of its quantum analog are displayed. As seen, although the Na atom is 
a relatively massive particle, at low temperatures the plateau is lower for the quantum 
case. This implies an initially faster decay arising from the strong influence of the quantum 
behavior at short time scales. It is therefore the real part of the intermediate scattering 
function what one should compare to the experiment rather than I2, as is usually done. 
Obviously, this effect will be less pronounced at high coverages. 
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FIG. 4: (Color online.) Classical intermediate scattering function (j45p for Na diffusion on a Cu(OOl) 
surface at 50 K and 100 K (black solid lines) and the real part of its quantum-mechanical analog 
(jSip (red/dark grey dashed lines). The surface coverage considered here is 0.028 and the diffusion 
along the azimut [100]. 

At very low temperatures, the Matsubara series should play a similar role like the extra 
term observed in the / and g functions like in the harmonic case. Thus, a new quantum 
correction should be added for very low temperatures. 

IV. CONCLUSIONS 



It is remarkable that, within the Markovian formalism presented here, the quantum inter- 
mediate scattering function, /i, is independent of the relative corrugation of the surface and. 
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at short times, also independent of the friction. At low surface temperatures, the Ii factor 
will be responsible for a higher contribution of the imaginary part of /, given by Eq. flSTj) . 
modifying substantially the response in the diffusion process. For relatively heavy particles 
and at very long times (diffusion time scales), operators in the I2 factor can be replaced 
by variables, since is very small. As far as we know, an exact quantum calculation for 
a corrugated surface is not possible and some approximations have to be invoked, e.g., the 
damped harmonic oscillator has been applied here in order to obtain close formulas. Of 
course, other different, alternative theoretical approaches can also be found in literature 



(see, for example, Refs. 
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531 ) but within the single adsorbate approximation. The theo- 



retical formalism that we propose here should also be very useful to avoid extrapolations 
at zero surface temperature when trying to extract information about the frustrated trans- 
lational mode. Diffusion experiments at low temperatures are very difficult to perform (or 
even unaffordable); for example, for the HeSE technique surface temperatures around 100 K 
can be attainable. However, the type of theoretical calculations needed in this formalism is 
easy to carry out and they would provide a simple manner to go to lower temperatures with 
quite reliable results, thus allowing to extract confident values of magnitudes such as friction 
coefficients and oscillation frequencies. By decreasing the surface temperature, quantum ef- 
fects are extended at higher values of time. Going from 100 K to 50 K, the time where 
the quantum dynamics is important increases from 0.07 ps to 0.15 ps. The standard prop- 
agation time for diffusion is greater than 400 ps. In our opinion, the limits of applicability 
of this quantum theory should be at coverages up to 12-16 % and around 50 K where the 
Matsubara series is still playing no role on the diffusion. Obviously, at lower surface temper- 
atures, the quantum character of the noise becomes more and more important. This type 
of conditions as well as the diffusion mediated by tunneling needs to be more investigated 
since new experimental results are being analyzed. 
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